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Abstract: I discuss the current status of the comparison between theoretical predictions 
j /* ! | and experimental data, relevant to the production of open charm and bottom quarks 

in photon-hadron and photon-photon collisions. I advocate the use of a formalism that 
■ matches fixed-order computations to resummed computations in order to make firm state- 

ments on heavy flavour production as described by perturbative QCD. 

do: tf 

CD 

Heavy flavour physics has been traditionally a challenging testing ground for the predic- o 
tions of perturbative QCD. Loosely speaking, we like to define a quark as heavy when its 
mass to is much larger than A QCD . This property entails the possibility of computing in 
perturbation theory the cross section for the production of an open heavy quark, which is 



1. Introduction 



not possible in the case of a light quark. It is customary, although not always accurate, 
to say that the mass of the quark sets the hard scale of the production process, and thus 



the relevant parameter to the perturbative expansion is a s (m). Furthermore, the condi- 
tion m 3> A QCD leads us to expect that the perturbative predictions are only marginally 
affected by power corrections and by contributions of non-perturbative origin, that we can 
not compute from first principles. 

There is no doubt that the top quark is a heavy quark; in fact, perturbative QCD does 
a fairly good job in describing its production mechanism. There is also a consensus on the 
fact that the bottom can be consistently treated as a heavy flavour. The case of the charm 
is borderline; it is difficult to list it together with u, d, and s quarks; on the other hand, we 
expect non-perturbative physics to have a non-negligible impact in this case, and we know 
that perturbative corrections are huge, since a s (1.5 GeV) ~ 0.3. 

In this paper, I shall not deal with top physics, and I shall concentrate on charm and 
bottom production in collider processes where at least one of the incoming particles is 
an on-shell photon. Thus, I shall not treat bottom production at hadronic colliders, nor 
charm and bottom production at fixed-target experiments. I shall only briefly remind the 

'Speaker. 
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reader that Tevatron data for the production of B mesons lie above the QCD predictions 
obtained with "default" parameter choices; however, QCD can predict the shape of p T 
spectrum, and does well for a few bb correlations. Also, the comparison between theory 
and data improves if 6-jets are considered. As far as fixed-target charm hadroproduction 
is concerned, it appears that perturbative QCD cannot reproduce the data, unless non- 
perturbative effects, such as p T -kick, are supplemented. Fortunately, these are moderate, 
and the agreement is, all in all, satisfactory. More details can be found in ref. [pj 



2. Perturbative computations 

Using the factorisation theorem, I write the cross section for the production of an inclusive 
open heavy quark Q in photon-hadron collisions as follows: 

do~Q = if } ® <»7i + E ® ff ] ® 

j ij 

where and /.^ are the parton densities in the photon and in the hadron respectively, 
and da 7 i and d&ij are the short-distance cross sections, computed in perturbation theory. 
At present, they have been computed to the next-to-leading order (NLO) accuracy in a s , 
which means a em a 2 s and a| respectively. One has to keep in mind that behaves 



asymptotically as a em /a s , and thus oIctq in eq. ( |2,ip is a series in a em a l s ; the truncation 
of such a series at i = 2 will be denoted in what follows as fixed-order (FO) NLO result. 
As is well known, the two terms in the RHS of eq. ( |2.lD are separately defined in terms of 
Feynman diagrams, but have no physical meaning, and they must always be summed in 
order to obtain sensible physical predictions. 

A similar factorisation formula holds in the case of photon-photon collisions: 

dag = da 77 + ^ /j 7) ® (dfryj + da h ) + ^ ® /j 7) ® c%. (2.2) 

3 ij 

It is clear that the second and third terms in eq. G2.2|) are analogous to the two terms in 



the RHS of eq. (2.1); on the other hand, the first term in the RHS of eq. (2.2) is peculiar of 
photon-photon collisions, and it corresponds to those events in which two pointlike photons 
initiate the hard scattering. Also in the case of eq. Q2.2] ) all the terms in the RHS must be 
summed in order to obtain measurable quantities. 

Although theoretically well defined, open heavy quark cross sections are not directly 
measurable; a description of the hadronization of the heavy quark into a heavy-flavoured 
hadron is necessary in order to compare theoretical predictions to data. This is done as 
prescribed by the factorisation theorem through the following equation: 

d 3 a H (k) /• , ,d 3 a Q (k) 



d 3 k 



f D NP (z) d aQ ^ 8 z {k - zk)d 3 k dz, (2.3) 
J d 3 k 



where H is the heavy-flavoured hadron with momentum k, and k is the momentum of 
the heavy quark. L>np(-z) is the non-perturbative fragmentation function, which is not 
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calculable but is universal; in what follows, I shall adopt Peterson form [||. In eq. (]2.3|), 
it is assumed that fragmentation scales the 3-momentum of the incoming quark. Different 
prescriptions are possible, but in all cases k remains parallel to k; the mass shell condition 
can be either k 2 = m 2 , or k 2 = m 2 H ] and, finally, none of these prescriptions is boost invari- 
ant. However, all possible prescriptions coincide in the large-p T limit. These ambiguities 
turn into uncertainties in the physical cross sections, which should be taken into proper 
account when comparing QCD predictions to data. A complete study on this issue will be 
presented elsewhere Q; here, I just state the fact that these uncertainties are negligible 
with respect to the uncertainties due to the dependence upon other input parameters, such 
as mass and scales. 



3. Charm production 

In this section, I compare FO NLO predictions to data for D* meson production at HERA 
(jp collisions) and at LEP (77 collisions). The relevant computer codes have been devel- 
oped in ref. @ (for jp collisions) and in ref. fH (for 77 collisions). I shall set m=1.5 GeV, 
and the renormalization scale equal to the transverse mass of the quark, m T = \Jp% + m 2 . 
The factorization scale will be set equal to m T in jp collisions, and equal to 2m T in 77 
collisions, since in the latter case smaller values of p T are probed. The parton densities in 
the proton are given by the CTEQ5M1 set. As far as the photon is concerned, I shall use 
the AFG set in the case of jp collisions, and the GRS set in the case of 77 collisions; AFG 

has been adopted since in photoproduction the formalism of ref. pfl is used, which requires 

(5) 

densities defined in the MS subtraction scheme. I shall set A Q q D = 226 MeV, as constrained 

by the CTEQ5M1 set; this value is almost identical to the central value of the PDG global 

fit. The probability of a c quark fragmenting into a D* meson is P c ^d* = 23.5%. The 

on-shell photons at HERA and LEP are emitted quasi-collinearly by the incoming leptons. 

Their spectrum can thus be described by the Weizsacker- Williams formula; here, I shall 

use the form of ref. @ . 

In order to define their photoproduction events, HI and ZEUS adopt different cuts on 

the fraction y of the electron momentum carried away by the photon, and on the virtuality 

Q 2 of the photon: 

ZEUS : 0.187 < y < 0.869, Q 2 < 1 GeV 2 , (3.1) 

HI : 0.29 < y < 0.62, Q 2 < 0.01 GeV 2 , (3.2) 

(in this paper, I shall only deal with data obtained by HI with the ETAG33 electron tag- 
ger). In fig. [I] I present the ratio of the data relevant to D*-meson p T over FO theoretical 
predictions. The spectra have been measured by HI [|| and ZEUS Q experiments in dif- 
ferent visible regions; apart from the differences already pointed out in eqs. (3.1) and (|3.2j), 



HI impose a cut on rapidity (\y\ < 1.5), whereas ZEUS impose a cut on pseudorapidity 
(\n\ < 1.5). For each set of data, I compute FO predictions for two values of the e param- 
eter entering the Peterson function (e = 0.02 and e = 0.036), in order to give an estimate 
of the uncertainties due to the choice of this parameter, as constrained by recent fits [10]. 
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Figure 1: Ratio of data over theory (FO NLO) for the p T spectrum of D* mesons, in the visible 
regions of the ZEUS and HI experiments. The weighted averages are also given. 



Shape-wise, the theoretical cross section in the visible region appears to be only moder- 
ately sensitive to the choice of the e parameter. The smaller e value gives a slightly better 
description of the data, although e = 0.036 is theoretically preferred when used in the 
context of a FO computation (see ref. flQ|| ). Regardless of the value of e, HI data appear to 
be in agreement with FO predictions, while ZEUS data display discrepancies. Taking the 
data at face value, the two data sets also indicate different p T spectrum shapes. It has to 
be stressed that ZEUS data have smaller statistical errors. Given the different conclusions 
on QCD predictions that can be drawn by looking at the results of the two experiments, 
it is impossible to issue a unique statement on the comparison between theory and data. 
Should this problem persist when more data will be available, it will be necessary to define 
the same visible cross section within the two experiments. 

The same pattern can be observed in the case of the rapidity /pseudorapidity spectra, 
presented in fig. ^ and fig. || for HI and ZEUS data respectively. HI data are in general 
statistically compatible with FO NLO predictions (obtained with e = 0.036), ZEUS data 
are not. Unfortunately, as in the case of the p T spectrum, the cuts imposed in order to 
define the distributions are different. ZEUS data seem to suggest a shape different from 
that predicted by QCD, which fails to describe the data especially in the positive-?] region. 
The last data point, however, is by far the one affected by the largest error. A similar 
trend can be possibly seen in HI data, for the cuts p T > 2.5 GeV and 3.5 < p T < 5 GeV, 
but in this case the discrepancy is not statistically significant. 

I now turn to the production of D* mesons in 77 collisions, which is measured by 
LEP experiments by applying a (possibly effective) anti-tag condition on the scattered 
electrons and positrons, < 9 max ] this condition can be translated in a suitable form of 
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Figure 2: Ratio of HI data over theory (FO NLO) for the y spectrum of D* mesons, in the case 
of different p T cuts. 
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Figure 3: Ratio of ZEUS data over theory (FO NLO) for the n spectrum of D* mesons, in the 
case of different p T cuts. 



the Weizsacker- Williams function 0. I shall compare here the FO NLO predictions to 
OPAL [11, 12] and L3 [13, 14 data; the two experiments have slightly different visible 
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Figure 4: Ratio of data over theory (FO NLO) for the p T spectrum of D* mesons, in the visible 
regions of the OPAL and L3 experiments. 
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Figure 5: As in hg. 0, for the n spectrum. 



regions: 

OPAL : 6 max = 0.033, 2 < p T < 15 GeV, \rj\ < 1.5, (3.3) 

L3 : e max = 0.030, 1 < p T < 12 GeV, |7?| < 1.4. (3.4) 

Also, the average center-of-mass energies relevant to the data of the two experiments are 
different: yfs e + e - = 193 GeV and ^/s e + e - = 198 GeV for OPAL and L3 respectively. 

The ratio of data over FO predictions is presented in fig. || and fig. [5] for p T and r\ 
spectra respectively. In this case, only e = 0.036 has been considered. By taking the data 
at face value there is a weak indication of a p T spectrum softer than the one predicted by 
QCD; as far as 77 spectrum is concerned, data seem to agree with QCD computations. All 
the data lie above the theoretical predictions; it can be shown that L3 data are within the 
band obtained by stretching the parameters entering the computation (mass and scales), 
while OPAL data lie just above the upper end of this band (see ref. ||). In spite of this 
small discrepancy, it is fair to say that QCD gives a reasonable description of the current 
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data. A firmer conclusions will be drawn when more data will be available; also in this 
case, the comparison of theory and data would benefit if similar visible regions were defined 
by the different experiments. 



4. Bottom production 

Bottom rates are much smaller than charm rates (about three order of magnitude), and 
it is painful for experiments to collect the statistics sufficient to perform a measurement. 
In spite of this, recent years have witnessed a great progress in this field, and quite a few 
experimental results are now available. Different experiments use different techniques, and 
the measured observables are rather inhomogeneous: they can be visible cross sections, 
the visible region being defined by means of cuts applied to the bottom quark variables or 
to the variables of the \x produced in the decay, or they can be total rates, obtained by 
extrapolating the visible rates to the whole phase space. For this reason, at present the 
only way to obtain a coherent picture is that of comparing the data to a given theory, in 
this case NLO QCD. This is what is done in fig. || where FO predictions are compared to 
HI 



, ZEUS |J, OPAL gj, and L3 @ data 
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Figure 6: Ratio of data over theory (FO NLO) for total bottom rates, as measured in photopro- 
duction, DIS, and 77 collisions. 

It is striking that, for all the measurements except one, the ratio data/theory exceeds 
3. Not only these values are much larger than those that we get at the Tevatron, they are 
also much larger than the corresponding results relevant to charm production, as measured 
by the same experiments. From the point of view of QCD, this is rather difficult to explain: 
as mentioned in the introduction, the mass of the quark sets the scale for the production 
process, and we would expect bottom cross sections to be predicted more accurately than 
charm cross sections. Thus, fig. ^ calls for an explanation; either a standard one (a better 
understanding of the fragmentation mechanism, or a more accurate description of semilep- 
tonic decays), or a more involved one (QCD processes may not be the only production 
mechanisms at work); and, of course, we need the statistics to be increased. 



5. Beyond fixed-order computations 

The fixed-order computations described in section work fine as long as all the mass scales 



International Europhysics Conference on HEP 



Stefano Frixione 



relevant to the problem are of the same order of magnitude as the hard scale that is used to 
compute a s . If this is not true, the coefficients of the expansion in a s can be numerically 
large, since they depend upon logQi/^j where Q\ and Q2 are two of the mass scales. In 
other words, the expansion parameter is not a s any longer, but rather a s logQ\/Q2- In 
charm physics, this situation is easily encountered, when the p T spectrum is measured; if 
p T 3> m, terms such as logp T /m grow large. Techniques exist that can take into account 
the dominant logarithmic terms to all orders in a s ; the resulting cross sections are denoted 
as "resummed" (RS, also improperly called "massless"). In heavy flavour photoproduction, 
currently the resummation has been performed to the next-to-leading logarithmic level; that 
is, all the terms of order a em a s (a s log p T /m) k and a em a 2 s {a s \ogp T /m) k , k = l,---,oo, 
are included in the RS cross section. 

As a rule of thumb, one would then compare data to FO predictions when p T ~ m, and 
to RS predictions when p T m. The problem is, the inequality p T 3> m cannot be turned 
into a quantitative statement. It is thus desirable to write the single-inclusive cross section 
in a form that is sensible in the whole p T range, that is, which interpolates between the 
FO result, relevant to the small- and intermediate-p T regions, and the RS result, relevant 
to the large-p T region. This is the aim of ref. [19] and ref. ||, relevant to hadro- and 



photoproduction respectively. The main results of these papers read as follows: 

FONLL = FO + (RS - FOMO) x G{m,p T ) , (5.1) 

where FONLL (for Fixed Order plus Next-to-Leading Logarithms) gives sensible predic- 
tions in the whole p T range, and FOMO is obtained from FO by letting to zero all the terms 
suppressed by powers of m/p T . The subtraction of FOMO from RS in eq. (|5.1| ) is necessary 
to avoid double counting, since some of the logarithms appearing in RS are already present 
in FO. To be more precise, FONLL has the following features: 

• All terms of order a em a s and a em a 2 s are included exactly, including mass effects; 

• All terms of order a em a s (a s logp T /m) k and a em a 2 s (a s logp T /m) k are included, with 
the possible exception of terms that are suppressed by powers of m/p T . 

Finally, the function G(m,p T ) is rather arbitrary, except that it must be a smooth function, 
and that it must approach one when m/p T — ► 0, up to terms suppressed by powers of m/p T . 
In what follows, we shall use 

G(m,p T )= 2 f\ 2 , (5.2) 
p T + c z m z 



with c = 5. The practical implementation of eq. ( |5.1| ) is rather involved, especially in the 
photoproduction case; all the details can be found in ref. Q. So far, no FONLL results are 
available for the case of 77 collisions. 



A detailed study of the phenomenological consequences of eq. (5.1) relevant to charm 
production at HERA will be presented elsewhere ||. Here, I shall only repeat the study 
performed in section |3|, presenting the ratio of HERA data for the p T spectrum over the 
FONLL results. For the latter, the Peterson parameter has been set equal to e = 0.02, 
consistently with the findings of ref. [1C]. The results are presented in fig. [7], which has to 
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be compared to fig. |l]. The plots are rather similar; the average of the values data/theory is 
only marginally larger in the case of FONLL computations, and this behaviour is basically 
driven by the large-p T points. At a first glance, this seems to be counterintuitive, since 
FONLL is expected to perform better than FO at large p T ; but actually, it means that, 
shape-wise, the p T spectrum predicted by FONLL agrees better with data with respect 
to that predicted by FO. However, it has to be stressed that this is not yet statistically 
significant. It also becomes clear that, although FONLL improves over FO, FO predictions, 
and not RS predictions, were the right choice so far to compare to experimental data. It 
does not make much sense to compare these data to pure RS predictions; the p T is simply 
not large enough. It is important to notice that any agreement between RS predictions 
and data in this p T range must be regarded as accidental, and QCD is actually not tested 
at all. 



6. Conclusions 

It does not come as a surprise that NLO QCD undershoots charm data, or at least it does 
so for a "default" choice of parameters. However, the agreement with the experimental 
results is reasonable. Of all the data considered here, those of ZEUS are the only ones that 
can not described even with an extreme choice of parameters. The p T spectrum measured 
by ZEUS is harder than that of QCD, and the rj spectrum grows faster than QCD predicts 
towards the positive 77's. LEP data only marginally favour a softer p T spectrum than NLO 
QCD. The increase of the statistics and the extension of the measurements to larger p T 's 
will shed further light on these issues. The whole p T range in photoproduction can now be 
consistently treated within a single formalism, reducing the ambiguities in the comparison 
between theory and data. The use of similar visible regions by the different experiments 
working at the same machine will also help in performing more stringent tests of the theory. 

The new results relevant to bottom production are quite puzzling; the rates are dra- 
matically larger than QCD predictions, the disagreement with theory being much worse 
than the corresponding one in charm production. This is in contradiction with the picture 
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of the hard production process at work in QCD. However, more work has to be done in 
this field, both by theorists and experimentalists, before firm conclusions can be reached. 
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